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I. INTRODUCTION 



One of the major goals of high-energy heavy-ion re- 
search is to explore properties of strongly interacting 
matter, particularly its phase structure [ij. A great body 
of experimental data has been already accumulated by 
now. Theoretical analysis of these data is still in progress. 
Explanation of transverse-mass spectra of hadrons (espe- 
cially kaons) in collisions of heavynuclei turned out to 
be one of the most difficult tasks [1, d, I3 ■ 

Experimental data on transverse-mass spectra of kaons 
produced in central Au-|-Au (sj and Pb-|-Pb 0] colli- 
sions reveal peculiar dependence on the incident energy. 
The inverse-slope parameter (so called effective temper- 
ature T) of these spectra at midrapidity increases with 
incident energy in the energy domain of BNL Alternat- 
ing Gradient Synchrotron (AGS) and then saturates at 
energies of CERN Super Proton Synchrotron (SPS). In 
Refs. [1, Q it was assumed that this saturation is as- 
sociated with the deconfinement phase transition. This 
assumption was indirectly confirmed by the fact that mi- 
croscopic transport models, based on hadronic degrees of 
freedom, failed to reproduce the observed behavior of the 
kaon inverse slope 0, Q • Hydrodynamic simulations of 
Ref. [3| succeeded to describe this behavior. However, in 
order to reproduce it these hydrodynamic simulations re- 
quired incident-energy dependence of the freeze-out tem- 
perature which almost repeated the shape of the cor- 
responding kaon effective temperature. This happened 
even in spite of using an equation of state (EoS) involving 
a phase transition into the quark-gluon plasma (QGP). 
This way, the puzzle of kaon effective temperatures was 
just translated into a puzzle of freeze-out temperatures. 
Moreover, results of Ref. Q imply that peculiar incident- 
energy dependence of the kaon effective temperature may 
be associated with dynamics of freeze-out. 
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In Ref. it was shown that dynamical description of 
freeze-out 11 1 , accepted in the model of 3-fluid dynamics 
(3FD) [l3, 13.M.[lq. naturally explains the incident en- 
ergy behavior of inverse-slope parameters of transverse- 
mass spectra observed in experiment. This freeze-out dy- 
namics, effectively resulting in a pattern similar to that 
of the dynamic liquid-gas transition, differs from conven- 
tionally used freeze-out schemes. In the brief letter [l3| 
only midrapidity inverse-slope parameters in central col- 
lisions were presented. In t he p resent paper we would like 
to extend analysis of Ref. [l^l by presenting transverse- 
mass spectra themselves at various impact parameters 
and rapidities and compare them with available experi- 
mental data in the range from AGS to SPS incident en- 
ergies. These spectra are computed with the same set of 
model parameters as that summarized in Ref. [l^ . In 
particular, the hadronic EoS [l6j with incompressibility 
K = 210 MeV is used. 



II. THE 3FD MODEL 

The 3FD model is designed for simulating heavy-ion 
collisions in the range from AGS to SPS energies. Un- 
like the conventional hydrodynamics, where local instan- 
taneous stopping of projectile and target matter is as- 
sumed, a specific feature of the dynamic 3-fluid descrip- 
tion is a finite stopping power resulting in a counter- 
streaming regime of leading baryon-rich matter. This 
counter-streaming is described in terms of two interacting 
baryon-rich fluids, initially associated with constituent 
nucleons of the projectile (p) and target (t) nuclei. In ad- 
dition, newly produced particles, populating the midra- 
pidity region, are associated with a baryon-free "fireball" 
(f) fiuid. 

We have started our simulations [l3i[l2j[l3| with a sim- 
ple hadronic EoS [iBl- The 3FD model turned out to be 
able to reasonably reproduce a large body of experimen- 
tal data [13, [H, ll3l in a wide energy range from AGS 
to SPS. This was done with a unique set of model pa- 
rameters summarized in Ref. [13 . Problems were met in 
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description of the transverse flow [TJ] . The directed flow 
required a softer EoS at top AGS and SPS energies (in 
particular, this desired softening may signal occurrence 
of the phase transition into the QGP). 

The transverse-mass spectra are most sensitive to the 
freeze-out parameters of the model. In the 3FD model 
the same freeze-out procedure is applied to both the 
thermal and chemical processes. Normalization of the 
meson spectra is directly related to the freeze-out tem- 
perature. Inverse slopes of the transverse-mass spectra 
represent a combined effect of the temperature and col- 
lective transverse flow of the hydrodynamical expansion. 
Had it been only the effect of thermal excitation, in- 
verse slopes for different liadronic species would approxi- 
mately equal. The collective transverse flow makes them 
different. These two effects partially compensate each 
other: the later freeze-out occurs, the lower temperature 
and the stronger collective flow arc. Nevertheless, in- 
verse slopes turn out to be sensitive to the instant of the 
freeze-out. Below we demonstrate that not only slopes 
but also normalizations of the pion transverse-mass spec- 
tra are well described by this uniflcd chemical-thermal 
freeze-out. Situation with kaons is more delicate, since 
we have to take into account the associative production 
of strangeness giving rise to a strangeness suppression 
factor. 

Freeze-out procedure adopted in the 3FD model was 
analyzed in detail in Ref. [ll|. This method of freeze- 
out can be called dynamical, since the freeze-out process 
here is integrated into fluid dynamics through hydrody- 
namic equations. The freeze-out front is not deflncd just 
"geometrically" on the condition of the freeze-out crite- 
rion met^ but rather is a subject the fluid evolution. It 
competes with the fluid flow and not always reaches the 
place where the freeze-out criterion is met. This kind of 
freeze-out is similar to the model of "continuous emis- 
sion" proposed in Ref. There the particle emission 
occurs from a surface layer of the mcan-free-path width. 
In our case the physical pattern is similar, only the mean 
free path is shrunk to zero. Wc would like to mention 
that recently the continuous emission model was further 
microscopically developed [H, [l^ . 

In particular, this dynamical freeze-out results in a pe- 
culiar incident-energy dependence of the actual freeze- 
out energy density averaged over space-time evolution of 
the colhsion, (£out), see Fig. [T] As seen, (eout) reveals 
fast rise at AGS energies and saturation at the SPS ener- 
gies. This happens in spite of the fact that our freeze-out 
condition involves only a single constant parameter — the 
"triger" freeze-out energy density Efrz = 0.4 GeV/fm'^ — 
which was taken the same for all incident energies with 
the exception of low incident energies, for which we 
used lower values: ehz{2A GcV) = 0.3 GeV/fm^ and 
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FIG. 1: (Color online) Actual average freeze-out energy den- 
sity in central Pb-|-Pb collisions as a function of invariant 
incident energy. 



efrz(l^ GcV) ~ 0.2 GeV/fm'^. In our previous paper [T^l 
we have performed only a rough analysis of this kind^. 
To find out the actual value of Eout, we analyzed results 
of actual simulations [Tlj . 

The "step-like" behavior of (Sout) is a consequence of 
the freeze-out dynamics, as it was demonstrated in Ref. 
PH . At low (AGS) incident energies, the energy density 
achieved at the border with vacuum, e", is lower than 
Efrz- Therefore, the surface freeze-out starts at lower en- 
ergy densities. It further proceeds at lower densities up to 
the global freeze-out because the freeze-out front moves 
not faster than with the speed of sound, like any pertur- 
bation in the hydrodynamics. Hence it cannot overcome 
the supersonic barrier and reach dense regions inside the 
expanding system. With the incident energy rise the en- 
ergy density achieved at the border with vacuum gradu- 
ally reaches the value of Efrz and then even overshoot it. 
If the overshoot happens, the system first expands with- 
out freeze-out. The freeze-out starts only when e" drops 
to the value of efrz. Then the surface freeze-out occurs 
really at the value « e^z and thus the actual freeze-out 
energy density saturates at the value (eout) ~ £frz/2; i.e. 
at the half fall from e* to zero. This freeze-out dynamics 
is quite stable with respect to numerics [Tlj . 

Fig. [2]clarifies this "step-like" behavior in terms of the 
average temperature, transverse velocity and baryon den- 
sity achieved at the freeze-out in central Au-|-Au (at AGS 
energies, impact parameter 5 = 2 fm) and Pb-fPb (at 



^ The freeze-out criterion demands that the energy density of the 
matter is lower than some value efj-^- 



^ This is why in the main text of Ref. [l3 we mentioned the 
value of approximately 0.2 GeV/fm^ for Eout and in the appendix 
explai: 
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explained how the freeze-out actually proceeded. In terms of Ref. 
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FIG. 2: (Color online) Average temperature (Tfrz, see scale 
at the l.h.s. axis), transverse velocity (/? = vt/c, see scale at 
the r.h.s. axis) and baryon density over the normal nuclear 
density (nf^^/'^o, see scale at the r.h.s. axis) at the freeze- 
out in central Au+Au (at AGS energies) and Pb+Pb (at SPS 
energies) collisions as a function of invariant incident energy. 
Arrows point to the proper axis scale. 



SPS energies, 6 = 2.5 fm) collisions. The freeze-out tem- 
perature Tfrz reveals a similar "step-like" behavior. At 
SPS energies the freeze-out temperatures in Fig. [5] arc 
noticeably lower than those deduced from hadron mul- 
tiplicities in the statistical model [lO, HH. The reason 
for this is as follows. Whereas the statistical model as- 
sumes a single uniform fireball, in the 3FD simulations at 
the late stage of the evolution the system effectively con- 
sists of several "fireballs" : two overlapping fireballs (one 
baryon-rich and one baryon-free) at lower SPS energies 
and three fireballs (two baryon-rich and one baryon-free) 
at top SPS energies At the top SPS energies these 
three fireballs turn out to be even spatially separated. 
Therefore, whereas high multiplicities of mesons and an- 
tibaryons are achieved by means of high temperatures 
in the statistical model, the 3FD model explains them 
by an additional contribution of the baryon-free fireball 
at a lower temperature. The freeze-out baryon density 
nfrz exhibits a maximum at incident energies of i^iab = 
1 OA-BOA GcV which arc well within range of the planned 
FAIR in GSI. This observation agrees with that deduced 
from the statistical model [1^, even baryon density val- 
ues in the maximum are similar to those presented in Ref. 

[13. 

This dynamics of the freeze-out is of key importance 
for explanation of the transverse-mass spectra. In partic- 
ular, the "step-like" behavior of inverse slopes (effective 
temperatures) in Fig. [3] is related to the similar "step- 
like" behavior of (eout) (see Fig. [1]) or freeze-out tem- 
perature Tfrz (Fig. (2). Of course, this relation is not 
direct, since the collective transverse velocity also con- 
tributes to effective temperatures. However, since the 
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FIG. 3: (Color online) Inverse-slope parameters of transverse- 
mass spectra of kaons, pions and protons at midrapidity pro- 
duced in central Au+Au and Pb-|-Pb collisions as a function 
of invariant incident energy. Solid lines correspond to purely 
exponential fit by Eq. (O with A = 0, while dashed lines, as 
well as experimental data at AGS energies, correspond to fits 
with A = — 1 for pions and with A = 1 for protons. Experi- 
mental data are from Kefs. [1, [10, Sill- 



collective transverse velocity at the freeze-out does not 
strongly change in the considered energy range (see Fig. 
[5]), effective temperatures closely follow the shape of the 
(£out) energy dependence. These inverse slopes T were 
deduced by fitting the calculated spectra by the formula 
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where my and y are the transverse mass and rapidity, 
respectively, and A is a parameter which is taken different 
in different experimental fits, see, e.g., 0, 
Numerical problems, discussed in Ref. 
us from simulations at RHIC energies. Already for the 
central Pb-|-Pb collision at the top SPS energy the code 
requires 7.5 GB of (RAM) memory. At the top RHIC en- 
ergy, required memory is three order of magnitude higher, 
which is unavailable in modern computers. 



III. PROTONS 

The 3FD model does not distinguish isotopic species. 
In particular, the proton and neutron are considered to 
be identical, i.e. the nucleon. Therefore, a spectrum 
of protons is just associated with the nucleon spectrum 
multiplied by the factor Z/A with Z and A being the 
charge and the mass number of colliding nuclei, respec- 
tively (in this paper we consider only identical colliding 
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m^-m [GeV] 

FIG. 4: (Color online) Transverse-mass spectra of protons in various rapidity bins (centered at \y — yc.m.\) from central Au+Au 
collisions at incident energies -Eiab = 2A, AA, QA, and SA GeV. 3FD results are presented for impact parameter b = 2 fm. 
Midrapidity spectra are shown unsealed, while every next data set and the corresponding curve (from top to bottom) are 
multiplied by additional factor 0.1. Data are from E895 Collaboration [5^: open squares for midrapidity, circles for negative 
{y ~ J/c.m.) and triangles for positive {y — yc.m.)- 
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m^-m [GeV] 

FIG. 5: (Color online) Transverse mass distributions of protons from Au(8A GeV)+Au collisions at various impact parameters 
(fe) and in Vcirioxis rctpidity bins (centered at y — yc.m. ,). The most backward rapidity spectra are shown unsealed, while every 
next data set and the corresponding curve (from top to bottom) are multiplied by the additional factor 0.1. The percentage 
shows the fraction of the total reaction cross section, corresponding to experimental selection of events. Experimental data are 
from E917 Collaboration 



nuclei) . This approximation is quite good at low incident 
energies, when the number of produced secondary parti- 
cles (mostly pions) is well smaller than A. The particle 
production tends to restore isotopic symmetry. Hence, 



the Z / A approximation becomes worse at high incident 
energies. 

Reproduction of available transverse-mass spectra of 
protons is presented in Figs. HHSl We consider here only 
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FIG. 6: (Color online) Inverse-slope parameters of transverse-mass spectra of protons produced in Au-|-Au collisions at incident 
energies -Eiab = 6A, 8A, 10. GeV at various centralities as a function of rapidity y — yc.m.- 3FD results are presented for 
impact parameters b — 2, 4, 6, 8, and 11 fm (from top row of panels to bottom one). The percentage indicates the centrality, 
i.e. the fraction of the total reaction cross section, corresponding to experimental selection of events. Experimental data are 
from E917 Collaboration [H. 
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FIG. 7: (Color online) Transverse-mass distributions of pro- 
tons at the midrapidity from central Pb-|-Pb collisions at in- 
cident energies ^lab = 158^1, 80A, AOA SOA and 20^4 GeV. 
3FD results are presented for impact parameter b = 2.5 fm. 
The spectrum for 158^4 GeV energy is shown unsealed, while 
every next data set and the corresponding curve (from top 
to bottom) are multiplied by the additional factor 0.1. NA49 
experimental data are taken from ^ {158 A GeV), ^ {80 A 
and 40.4 GeV) and [3 {30 A and 20 A GeV). 



collisions of heavy nuclei, since they offer favorable con- 
ditions for application of the hydrodynamics. With few 
exceptions, overall reproduction of these spectra is quite 
good in terms of both normalization and slope. 

The first exception is the spectrum at the lowest con- 
sidered energy of 2A GeV, see Fig. S] The calculated 
spectrum turns out to be steeper than the experimen- 
tal one in all rapidity bins. This spectrum is calculated 
with our default choice of the freeze-out energy density 
efrz(2A GeV) = 0.3 GeV/fm'^ accepted in our previous 
papers [13, [13, IS] • Variation of Sfrz in the range from 0.1 
to 0.6 GeV/fm-^ does not noticeably affects the slope of 
the spectrum. This fact agrees with the above discussed 
dynamics of the freeze-out at low energies. Indeed, the 
surface freeze-out occurs at lower energy density which 
is the same independently of £frz. Only global frceze- 
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FIG. 8: (Color online) Transverse-mass distributions of pro- 
tons in the rapidity range 2.4 < y < 2.8 from Pb-|-Pb colli- 
sions at incident energy -Ei^b = 158 A GeV and various cen- 
tralities, i.e., impact parameters b. The spectrum for central 
collisions is shown unsealed, while every next data set and the 
corresponding curve (from top to bottom) are multiplied by 
the additional factor 0.5. NA49 experimental data are taken 
from Ref. [S^]. 3FD results are presented for rapidity y — 2.6 
and impact parameters estimated in Ref. [2^ . 



out of the system residue in the very end of the system 
disintegration produces slight sensitivity to efe. This 
problem at the energy of 2A GeV indicates that the ac- 
cepted model of the freeze-out is poorly applicable at low 
incident energies. 

The 3FD model does not exhibit deviation from the 
exponential fall-off at (tot — to) < 0.2 GeV, observed 
in the experiment, as it is most clearly seen in Figs. [3] 
and [8] The same feature of the hydrodynamic calcula- 
tion was earlier reported in Ref. [29| . As it was shown 
[2^ . a post-hydro kinetic evolution (afterburner) is re- 
quired to produce the observable two-slope structure of 
the TOr-spectra. In our model such a post-hydro evolu- 
tion is absent. It is worthwhile to note that the above 
problem is not an inevitable feature of any hydrodynamic 
calculation. For instance, the deviation from the expo- 
nential fall-off was reproduced in calculations by Kolb et 
al. [13. 

Traditionally, this deviation from the exponential fall- 
off at low (tot — rn) is associated with collective trans- 
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FIG. 9: (Color online) Transverse- mass spectra of positive and negative pions in various rapidity bins (centered at \y — yc.m.\) 
for central Au+Au collisions at incident energies -Eiab = 2A, 4A, QA, and 8A GeV. 3FD results (for (tt"*" + tt" + 7r~)/3) 
are presented for impact parameter 6 = 2 fm. Midrapidity spectra are shown unsealed, while every next data set and the 
corresponding curve (from top to bottom) are multiplied by the additional factor 0.1. Data are from the E895 Collaboration 
[3^ : squares for midrapidity n~ , triangles for positive-rapidity 7r~, circles for negative-rapidity tt~ , diamonds for midrapidity 
and positive-rapidity vr"'", inverted triangles for negative-rapidity tt'^ . 
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FIG. 10: (Color online) Transverse-mass spectra of positive 
pions in various rapidity bins from central Au+Au collisions 
at incident energy -Ei^b = 10. GeV. 3FD results (for (tt"*" + 
7r°-|-7r~)/3) are presented for impact parameter b = 2 fm. The 
most backward rapidity spectrum is shown unsealed, while 
every next data set and the corresponding curve (from top to 
bottom) are multiplied by the additional factor 0.1. Data are 
from the E802 Collaboration [U. 



verse expansion of the system [2^, |3l| . In Ref. [1^ it is 
described how the post-hydro cascade modify this radial 
flow. Apparently the final-stage (post-hydro) Coulomb 
interaction also distinctly affects the flow, accelerating 
or decelerating it depending on the electric charge of the 
species. It is clear from the difference between spectra 
of positive and negative pions, see Fig. [HI As seen from 
Fig. [21 negative pions arc decelerated (their low-my spec- 
trum is enhanced) while positive pions arc accelerated 
(their Iow-tot spectrum is suppressed). The same mech- 
anism should also contribute to suppression of Iow-tot 
spectrum of protons as compared to pure hydrodynamic 
calculation. With rapidity moving off the midrapidity, 
mr-spectra of protons reveal weaker suppression of their 
low-mT parts, see Fig. 21 This fact complies with the 
Coulomb mechanism. Indeed, faster particles (in the 
frame, where the fireball generating the field is at rest) 



spend shorter time in the Coulomb field and, hence, ac- 
quire smaller additional momentum. 

Thus, keeping in mind that reproduction of the data 
in the Iow-tot region can be cured by application of a 
post-hydro dynamics, we see that, in general, the 3FD 
model reasonably reproduces proton m^-spectra beyond 
extreme kinematic regions — rapidities being to far from 
the midrapidity and too high ttit (see Fig. |4]) — in not too 
peripheral collisions of heavy nuclei (see Fig. [5]). In the 
above marginal regions applicability of hydrodynamics is 
questionable. However, even in very peripheral collisions 
agreement with data can be unexpectedly good, see Fig. 

m 

Inverse-slope parameters [cf. Eq. ([T]) with A = 1] , dis- 
played in Fig. [HI summarize our results at AGS energies. 
The most pronounced disagreement with experimental 
temperatures occurs in midrapidity region, especially at 
10.8^ GeV incident energy. This problem is directly re- 
lated to the above-discussed poor reproduction of low- 
rriT suppression of the spectra. Apart from this problem, 
effective temperatures arc well reproduced. 



IV. PIONS 

Comparison of calculated transverse-mass spectra of 
pions with available data is presented in Figs. [9]- 
[T2l Since the 3FD model does not distinguish isotopic 
species, we assume that numbers (and spectra) of tt"'", 
TT~ and tt" are identical and equal 1/3 of the total num- 
ber (spectrum) of pions. In view of this approximation, if 
calculated spectra are situated between measured spec- 
tra of positive and negative pions, that would mean a 
good agreement of our calculations with data. Contrary 
to protons, this approximation becomes better applicable 
at high incident energies, when abundant production of 
secondary particles results in partial restoration of iso- 
topic symmetry of the system. As it is illustrated in Fig. 
[21 spectra of positive and negative pions indeed become 
more and more similar with the incident energy rise. 

Similarly to the proton greement with data is 

the worst at the lowest considered incident energy of 2A 
GeV, see Fig. [HI There our calculated spectra closely 
follow experimental spectra for positive pions, while they 
should lie in between tt"*" and tt~ spectra. Above this en- 
ergy, overall reproduction of pion spectra is quite good 
both in normalization and slope. It is even better than 
that for protons. The calculated pion spectra also repro- 
duce the rapidity (Figs. [5] and [TH)) and incident-energy 
(Figs. [11] and [T^ dependence of the data. Note that 
this agreement is achieved with the unified chemical- 
thermal freeze-out. Of course, two different freeze-out 
criteria for chemical and thermal processes could slightly 
improve this agreement with the experiment. However, 
that would be too high price for so tiny improvement. 



10 




0.5 1 

Ytij-m [GeV] 

FIG. 11: (Color online) Transverse-mass spectra of positive 
pions at midrapidity in central Au+Au collisions at incident 
energies Ejab = 2A, 4A, 6A, 8A, and 10.7 A GeV. 3FD results 
(for (tt^ + tt" + 7r~)/3) are presented for impact parameter 
6 = 2 fm. Spectrum for 10.7^ GeV is shown unsealed, while 
every next data set and the corresponding curve (from top to 
bottom) are multiplied by the additional factor 0.1. Data are 
from the E866 Collaboration 0] . 
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FIG. 12: (Color online) Transverse-mass spectra of negative 
pions at midrapidity in central Pb-|-Pb collisions at various in- 
cident energies E^^^, = 20A, 30A, iOA, 80A, and 1584 GeV. 
3FD results (for (tt"*" -|- tt" -I- 7r~)/3) are presented for impact 
parameter b = 2.5 fm. The spectrum for 158A GeV energy 
is shown unsealed, while every next data set and the corre- 
sponding curve (from top to bottom) are multiplied by the 
additional factor 0.1. Data are from the NA49 Collaboration 
0- 



V. KAONS 

Upon demonstrating that calculated proton and pion 
transverse-mass spectra are in reasonable agreement with 
available data wc proceed to kaon spectra. The "step- 
like" incident energy dependence of the the inverse-slope 
parameter (T, cf. Eq. ([1])) of these spectra was inter- 
preted as a signal of the deconfinement phase transition 
1^, Q. Moreover, microscopic transport models, based 
on hadronic degrees of freedom, failed to reproduce the 
observed behavior of the kaon inverse slope 0, @| • 

Comparison of calculated of transverse-mass spectra of 
kaons with available data in collisions of heavy nuclei is 
presented in Figs. I13II16I Overestimation of kaon nor- 
malization at low incident energies (i?iab < 6A GeV) is 
not surprising, see Fig. [T31 This is not a result of the uni- 
fied chemical-thermal freeze-out. The reason is that we 
use the EoS based on grand canonical ensemble, therefore 
production of rare particles should be overestimated at 
low energies. This normalization can be easily corrected 



by introducing a suppression factor 7s, taking into ac- 
count associative production of strangeness. Apart from 
this normalization, overall reproduction of kaon spectra 
is quite satisfactorily. Separate treatment of the chemi- 
cal and thermal freeze-outs could certainly improve the 
normalization of kaon spectra at SPS energies. However, 
to do this we would have to introduce different freeze- 
out criteria for K'^ and K~ . This is distinctly seen, e.g., 
at the example of and spectra at i^iab — 80A 
GeV, see Figs. [15] and [1^1 The K'^ spectrum requires 
somewhat earlier freeze-out in order to increase the nor- 
malization. While the K~ spectrum demands for later 
freeze-out to reduce the normalization. In fact, there are 
physical arguments in favor of this sequence of K'^ and 
K~ freeze-outs. We prefer a cruder description of observ- 
ables to introduction of additional fitting parameters. 

Dependence on the rapidity is reproduced, see Fig. 1131 
as well as evolution of slopes with the incident energy 
variation, see Figs. [T5HTB1 Fig. [3] demonstrates that 
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m^-m [GeV] 

FIG. 13: (Color online) Transverse-mass spectra of positive 
kaons in various rapidity bins in central Au+Au collisions at 
incident energy i?iab = 10. 6^4 GeV. 3FD results are presented 
for impact parameter 6=2 fm. The most backward rapidity 
spectrum is shown unsealed, while every next data set and 
the corresponding curve (from top to bottom) are multiplied 
by the additional factor 0.1. Data are from the E802 Collab- 
oration [33l |. 




0.5 1 

nij-m [GeV] 



FIG. 14: (Color online) Transverse-mass spectra of positive 
kaons at midrapidity in central Au-|-Au collisions at incident 
energies i^jab = 2A, 4A, 6A, 8A, and 10.7 A GeV. 3FD resuhs 
are presented for impact parameter b — 2 fm. Spectrum for 
10.7j4 GeV is presented unsealed, while every next data set 
and the corresponding curve (from top to bottom) are mul- 
tiplied by the additional factor 0.1. Data are from the E866 
Collaboration Q|. 



VI. SUMMARY 



inverse slopes of kaon spectra are indeed reasonably re- 
produced. This inverse slopes were deduced by fitting 
the calculated spectra by formula H]) with A = (purely 
exponential fit). Moreover, the pion and proton effec- 
tive temperatures also reveal saturation at SPS energies, 
if they are deduced from the same purely exponential fit 
with A = 0. Though the purely exponential fit with A = 
does not always provide the best fit of the spectra, it al- 
lows a systematic way of comparing spectra at different 
incident energies. In order to comply with experimental 
fits at AGS energies (and hence with displayed experi- 
mental points displayed by squares), we also present re- 
sults of fits with A = — 1 for pions and with A = 1 for 
protons. The most pronounced disagreement with ex- 
perimental temperatures takes place for protons, while 
reproduction of the proton spectra themselves (see Fig. 
[7]) docs not look bad. 



Transverse-mass spectra of protons, pions and kaons 
in wide ranges of incident energies (from AGS to SPS), 
rapidity bins and centralities have been analyzed within 
the 3FD model. These spectra were computed with the 
same set of model parameters as that summarized in Ref. 
[T3 |. In particular, the hadronic EoS [l^ with incom- 
pressibility K = 210 MeV was used. We considered here 
only collisions of heavy nuclei, since they offer favorable 
conditions for application of the hydrodynamics. It was 
demonstrated that with few exceptions the 3FD model 
reasonably reproduces these spectra in the range of in- 
cident energies 4A GeV < -Biab < 160A GeV. The ex- 
ceptions concern extreme kinematic regions — rapidities 
being to far from midrapidity and too high niT — and 
too peripheral collisions of heavy nuclei. In the above 
marginal regions applicability of the hydrodynamics is 
questionable. However, even in very peripheral collisions 
agreement with data can be unexpectedly good. 
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nij-m [GeV] 



FIG. 15: (Color online) Transverse-mass spectra of positive 
kaons at midrapidity in central Pb+Pb collisions at various 
incident energies Bj^b = 20A, 30A, 40A, 80A, and 158A GeV. 
3FD results are presented for impact parameter b = 2.5 fm. 
The spectrum for 158A GeV energy is shown unsealed, while 
every next data set and the corresponding curve (from top to 
bottom) are multiplied by the additional factor 0.1. Data are 
from the NA49 Collaboration 0]. 



The 3FD model reproduces inverse slopes of rriT- 
spectra, in particular, the "step- like" dependence of kaon 
inverse slopes on the incident energy. This hydrodynamic 
explanation of the transverse-mass spectra and "step- 
like" behavior of effective temperatures implies that at 
considered incident energies a heavy nuclear system re- 
ally reveals a hydrodynamic motion during its expansion. 

In fact, microscopic transport and hydrodynamic mod- 
els do not contradict each other. The fact that they all 
successfully describe many observales in collisions heavy 
nuclei in spite of big differences between these models 
suggests that detailed information on cross sections is 
unimportant for these reactions. The collision frequency 
is already high enough in order to make hydrodynam- 
ics to be applicable. This conjecture has already been 
put forward in Ref. 34 1. In particular, transport 0, 0] 



and hydrodynamic models well describe pion transverse- 
mass spectra of pions in wide range of incident energies. 
Therefore, the failure of microscopic transport models 
0, Q to reproduce the observed behavior of the kaon 
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FIG. 16: (Color online) The same as in Fig. 1 151 but for nega- 
tive kaons. 



inverse slopes may be interpreted as a signature that 
kaon and antikaon interaction cross sections in micro- 
scopic models are not big enough for kaons to be cap- 
tured by matter in a common flow. In Ref. [2l a kind 
of such cross-section enhancement (Cronin effect) was 
tested against the experimental data. The resulting ef- 
fect of that enhancement was found insufficient to ex- 
plain the observed discrepancy. Another possibility is 
that multi-body collisions are important in the transport. 
This was checked within the Giesscn Boltzmann-Uehling- 
Uhlenbeck (GiBUU) model [35|, where three-body inter- 
actions were included in simulations. It was found that 
the three-body collisions indeed result in good reproduc- 
tion of all transverse-mass spectra [35j . 

Returning to the question if the considered " step- like" 
behavior of effective temperatures is a signal of a phase 
transition into the QGP, we have to admit that this is 
not quite clear as yet. It depends on the nature of the 
freeze-out parameter ehz = 0.4 GeV/fm'^ which should 
be further clarified. The EoS is not of prime importance 
for this behavior. The only constraint on the EoS is that 
it should be in some way reasonable. Moreover, our pre- 
liminary results indicate that a completely different (from 
that used in this work) EoS [11] with 1st order phase 
transition to the QGP still reasonably reproduces this 
"step-like" behavior even in spite of that it fails to de- 
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scribe a large body of other data. This happens because 
the same freeze-out pattern is accepted there. 

This kind of the freeze-out has its imphcations for 
observation of fluctuations. A number of observables 
were suggested, which are based on various fluctuations: 
transverse- momentum fluctuations [sj, electric-charge 
fluctuations (ssl . [s^ , baryon-charge-strangcncss correla- 
tions t40j], fluctuations of the K/ir ratio [41j; for the re- 
cent review see Ref. [l^. These observables distinguish 
between the hadronic and quark-gluon phases. However, 
available experimental data on fluctuations: on fluc- 
tuations [43, [H, [isl. 46l and positive-to-negative charge- 
ratio fluctuations 47l. l48l| , at SPS energies well agree with 
the hadronic predictions. Only the K/t: fluctuations [i^ 
reveal a nontrivial behaviour with the incident energy. 

The fact that the freeze-out takes place at rather low 
values of energy density and temperature (see Figs. [1] 
and 12) explains why we observe only hadronic fluctua- 
tions at the freeze-out, even if the compressed matter 
was in the quark-gluon phase. All the quark-gluon fluc- 
tuations get dissolved during the long path of the sys- 
tem from hadronization to freeze-out. Fluctuations of 
the K /-K ratio can be an exception from this rule, since 
this signal is the most robust. A meson can disap- 
pear only if it meets a particle of the opposite strangeness 
(A, S or K~). If the strange system is relatively di- 
lute at the hadronic stage of the expansion, as it is 
the case at lower SPS energies, the probability to meet 
a particle of the opposite strangeness is relatively low. 
Therefore, the signal from the quark-gluon phase par- 



tially survives. At higher energies (top SPS energies) 
the strange system becomes already dense enough to de- 
stroy this signal. This mechanism can explain a non- 
trivial energy dependence of the Kh: fluctuations, ob- 
served by the NA49 Collaboration [49|. It implies vio- 
lation of the chemical equilibrium for strange particles 
at the late stage of the expansion. Therefore, it is be- 
yond the frame of the present version of the 3FD model, 
where meson abundance immediately follow changes 
of densities of the matter. If this explanation is true, 
fluctuations of the ratio ii'+/7r+ = Ni(+/N^+ should re- 
veal even stronger energy dependence than those of the 
ratio K/tt = {Nk+ + Nk~)/{N^+ + N^-), where Na is 
the observed yield of a-particles. The same mechanism 
can be behind the horn effect in the excitation function 
of the K+/n+ ratio 
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